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We estimate fermion loop corrections to the two-point correlation function of primordial tensor 
perturbations in a slow-roll inflationary background. We particularly compute an explicit term of 
one-loop correction from a massless fermion, and then extend to the complete Interaction Hamil- 
tonian. After that, we study one-loop corrections contributed by a massive fermion to primordial 
tensor fluctuations. The loop correction arisen from a massless fermion field contains logarithms and 
thus may constrain the validity of perturbation theory in inflationary cosmology, but the situation 
could be relaxed once the fermion's mass is taken into account. Another one-loop diagram for a 
massive fermion which involves one vertex is constrained by a UV cutoff as expected by quantum 
field theory. Our result shows that loop corrections of a fermion field have the same sign as those 
of a scalar field, and thus implies that the inclusion of fermion loop corrections may not help to 
alleviate the issue of IR divergence in inflationary cosmology. 

^ PACS numbers: 98.80.-k, 98.80.Cq 

I. INTRODUCTION 

Inflation is considered as the most successful model of describing physics of very early universe, which has explained 
conceptual issues of Big Bang cosmology [lHH (see 0-Q for early works). Among these remarkable achievements, 
inflation has predicted a nearly scale-invariant primordialpower spectrum which was later verified in high precision 
by Cosmic Microwave Background (CMB) observations Q. The success of inflation is mainly based on a scries of 
assumptions including an enough long period of quasi-exponential expansion and the applicability of perturbation 
theory during this phase. Therefore, an important question which ought to be understood is whether there exists a 
bound beyond which the perturbation theory would break down in inflationary cosmology. 

Generally, the primordial perturbations during inflation were originated from quantum fluctuations with their phys- 
ical wavelengths being much shorter than the Hubble radius for which the modes are in ultraviolet (UV) regime. Along 
with the quasi-exponential expansion, the wavelengths of these modes would be stretched outside the Hubble radius 
■^j- ' for which the modes are in infrared (IR) regime. Through this mechanism quantum fluctuations are able to become 
classical ones and finally can be observed in today's CMB experiments (see |8| for an overview of cosmological pertur- 
bation theory). Therefore, the evolution of metric perturbation can be calculated at the moment of Hubble crossing 
\ until the end of inflation. This mechanism applies to both the scalar type and tensor type of metric perturbation 
1— I ■ during inflation. Especially, for tensor fluctuations the amplitude of one degree of freedom is determined in order of 
H/M p , and this result is not sensitive to the details of inflation models. 

Since the observable modes of inflationary fluctuations are in IR regime, the reliability of perturbation theory 
expects that these modes may behave well and two-point correlation functions may be finite even when higher order 
corrections are included. However, as was noticed by a pioneer work [j| long time ago, perturbation theory beyond 
leading order had troublesome IR behavior in a de Sitter space. Another issue of IR effect which may affect late time 
evolution of the universe due to the accumulation of long- wavelength inflationary fluctuations was studied in (lol - 13 j . 
In the past ten years, following Maldacena's calculation of three-point functions in a slow- roll inflation model 14|, 
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there have been significant theoretical developments on refining our understanding of quantum field description of 
inflationary fluctuations. A remarkable work systematically studying quantum effect to arbitrary order in cosmological 
perturbations was recently formulated by Weinberg [l5| (see [l6| for further studies of the time dependence of the 
correlator.). With an example of massless scalar loop corrections, it was observed that the momentum dependence 
of two-point correlation function of inflationary fluctuations evolves as a logarithmic function in IR limit [15| . This 
result implies that perturbation theory may break down since of the contribution of IR fluctuations to loop integrals 
fl7l - l20j . There have been a lot of studies on the effects of IR divergences on scalar type curvature perturbations in 
inflation (see [2l[ for a recent review; and see [22l - l38j for detailed studies of various inflationary loop corrections; and 
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also see [39l444| for fully dimensionally regulated and renormalized computations in inflation.). Note that, the IR 
effect of primordial inflationary perturbations in the presence of only adiabatic modes could be very limit [4514481 ). 
Therefore, the study of IR issue is much more interesting when entropy modes are taken into account, namely, a 
second scalar or a fermionic field minimally coupled to inflationary background. 

In this paper, we estimate the IR issue arising from loop corrections of a fermionic field to the two-point correlation 
function of primordial tensor fluctuations in slow-roll inflationary cosmology. We first consider loop contribution from 
a massless fermion. As an exercise, we calculate one term of the Interaction Hamiltonian and find that it contributes 
to one loop involving two vertices and gives rise to a logarithmic correction to the two-point correlation function of 
tensor fluctuations. After extending to the complete form of the Interaction Hamiltonian, this logarithmic correction 
remains while only the coefficient in front of it is improved. However, once we consider the fermion's mass is not zero, 
we find the logarithmic function would be replace by a cosh function and thus the IR issue could be much relaxed. 
In addition, for a massive fermion there exists another one-loop diagram which only involves one vertex, but this 
diagram only contributes to the UV regime and thus is not of enough observable interest. Our main result shows 
that the loop corrections contributed by a fermion field have the same sign as those contributed by a scalar field, and 
thus it implies that the inclusion of fermion loop corrections cannot be used to alleviate the issue of IR divergence in 
inflationary cosmology. 

The present paper is organized as follows. Section II gives a brief introduction to the background dynamics of 
a fermion field in inflationary cosmology. In Section III, we first simply review the canonical quantization of linear 
fluctuation modes in cosmological background, and then introduce the in-in formalism as a setup for computing loop 
diagrams. After that, we expand the inflationary Lagrangian involving a fermionic field to the third and forth order 
and thus derive the interaction Hamiltonian for a fermionic field coupled to tensor fluctuations straightforwardly. 
Making use of the in-in formalism and the interaction Hamiltonian, we in Section IV estimate loop corrections of 
both a massless and a massive fermionic field to the two-point correlation function of primordial tensor fluctuations 
up to one-loop order. Finally, we conclude with a discussion in Section V. Note that, we use the natural units 
8ttG = H = c = 1 and all parameters are normalized by M p = 1 / V 8ttG in this paper. 



II. BACKGROUND DYNAMICS 



To begin with, we simply review the background dynamics of a sp inor field which is minimally coupled to Einstein's 
gravity (see Refs. [4!|[5(| for detailed introduction and see [HTl l52| for recent phenomenological study in cosmology). 
Following the general covariance principle, a connection between the metric g^ v and the vierbein is given by 

9^ e a e b = Vab , (1) 

where denotes the vierbein, g^ is the space-time metric, and i] a b is the Minkowski metric with r] a \, = 
diag(I, —I, —I, —I). Note that the Latin indices represents the local inertial frame and the Greek indices represents 
the space-time frame. 

We choose the Dirac-Pauli representation as 

M\1X -<•-(! ?). m 



— 1 y ' ' V _<r * , 

where is Pauli matrices. One can see that the 4x47° satisfy the Clifford algebra {7 a ,7 b } = 2rj a b. The 7° and 
provide the definition of a new set of Gamma matrices 

I* = e£ 7 a , (3) 



which satisfy the algebra {r M , T 1 '} = 2g^ v . The generators of the Spinor representation of the Lorentz group can be 
written as S afc = j [7 , j b ] . So the covariant derivative are given by 

= (0„ + , (4) 

= d$ - , (5) 

where the Dirac adjoint is defined as 'tp = The 4x4 matrix f2 M = ^ui lia b'S, ab is the spin connection, where 

Wfj, a b = e^V^e^b are Ricci spin coefficients. 

By the aid of the above algebra we can write down the following Dirac action in a curved space-time background 

d 4 x e [l(^r"D M ^ - D^T^) - m^} . (6) 
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Here, e is the determinant of the vierbein e° and m stands for the mass of the spinor field ip. 
Varying the action with respect to the vierbein e^, we obtain the energy-momentum-tensor, 

1 _ _ _ _ 

= -[ipTvD^ip + ipV^D^tp - DpipTvip - D^T^ip] - g^Cp . (7) 

On the other hand, varying the action with respect to the field ip, ip respectively yields the following equations of 
motion, 

iT^D^ -mip = , (8) 
iD^V +m>tp = . (9) 

We deal with the homogeneous and isotropic FRW metric, 

ds 2 = dt 2 -a 2 (t)dx 2 . (10) 



Correspondingly, the vierbein are given by 



a 



Assuming the spinor field is space- independent, the equation of motion reads i"f°(ip + ^Htp) — mip = 0, where a dot 
denotes a derivative with respect to the cosmic time and H is the Hubble parameter. Taking a further derivative, we 
can obtain: 

N 

W = ^ , (12) 

where N is a positive time-independent constant. 

From the result (|12[) . we can find the background dynamics of a massive spinor field behaves like the cold dark 
matter with its energy density goes as a~ 3 . If the universe is in inflationary phase, the energy density of a spinor field 
will be diluted immediately and cannot contribute to the background evolution. Therefore, in the quantum treatment 
of cosmological perturbations in later context, we can safely choose (ip) = as the vacuum state for the fermion field 
during inflation. 

III. SETUP OF PERTURBATION ANALYSIS 

Inflationary cosmology usually assumed a homogeneous scalar field with a non-zero expectation value varying slowly 
along its potential. This background dynamics is characterized by a series of slow-roll parameter, among which the 
most significant parameter is given by e = —^jz- In this scenario, we naturally have cosmological perturbations of both 
scalar and tensor types. The interaction between the scalar perturbation and a fermion was studied comprehensively 
in [24l |. In the current paper, we particularly focus on the interaction between tensor fluctuation and a fermion field. 

A. Canonical quantization of linear fluctuations 

Taking into account the spatial dependence of the fermion field, we can write down the equations of motion for the 
fermion and tensor fluctuations as follows, 

h°(tp + -Hib)-i^— tb -mtp = , (13) 
2 a 

d 2 

h l3 +'mh l3 - -\hij = . (14) 
Moreover, we make the Fourier transformations of these fields, which are given by 

^(x,t) - J d 3 p^ [e lp - x a PtS X PiS (t) + e-^*(3l s W PtS (t)] , (15) 
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and 

%(x,i) = / d 3 g ]T [e^e lJ (q,X)a(q,X)h q (t) + e"'^ (g, A)a* (q, X)h* q (t)] , (16) 

J A 

where s = ±i is the spin of the fermion, A = ±2 stands for the helicity of tensor fluctuation, and eij(q,X) is the 
polarization of tensor mode with the following normalization [53j : 

^ e i j{q,X)e* H {q, A) = 5 ifc <5^ + S a S jk - <%4; + 8 i:j q k qi + 

A=±2 

-SikQjqi - S.aq j q k - S jk q. t qi - 5jiq t q k + qiqjQkqi , (17) 

where q is the unit vector along the q direction. 

The annihilation operators a P , s , /?p.s, ct{q, A) satisfy the following commutation relations 

{a(p, a), at(p', a')} - {/3(p, *), /3 f (p', *')} - ^8 3 (p - p') , (18) 
{a(p, «), o(p', s')} = {0(P, *), j8(p', ^')} - , (19) 

and 

[«(q,A), a *(q',A')l - <W 3 (q-q') , (20) 
[a(q,A),a(q',A')] = , (21) 

respectively. 

Additionally, the Fourier modes X PtS (t), W PtS (t) and h q (t) satisfy the equations of motion for the fermion and tensor 
fluctuation in cosmological background as shown in Eqs. p^|) and (fT4")) respectively. At quadratic order, each Fourier 
mode evolves independently. Therefore, we can solve the equations of motion for those fluctuations in comoving time 
coordinate r = J dt/a, and obtain 



Mr)-^^(l + «rK*. (23) 

for the fermion and tensor mode, respectively. Note that, in above solutions we have made use of Bunch-Davics vacuum 
as initial condition to determine the coefficient of the Hankel function, and the convention of primordial tensor modes 
is the same as that adopted in (54[. Also we notice that, when the fermion is massless (m=0), the solutions to the 
fermion modes are independent of the helicity and take the form of ^(^a) 3 / 2 , which is asymptotically the plane wave 
function. 

B. Interactions 

After obtaining solutions of linear fluctuations, we then study the interaction terms for loop calculation. Since 
the Hamiltonian that governs primordial fluctuations is time dependent in our case, we need to make use of in-in 
formalism. Following Weinberg's formula (l5j . we can calculate the in-in correlation function to arbitrary order via 
the following expression, 

(Q(t)) = Yi N dt N dt N -i... dhdH!^),^!^),...^!^)^ 1 ^)}...}}) , (24) 



^=0 
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where Hi is the interaction Hamiltonian as will be discussed later. In this formula, the LHS is working in the interactive 
vacuum and the expectation value on the RHS is calculated in the effectively free field vacuum by annihilation operators 
working on cosmological background. Moreover, the index "N" in Eq. (|24[) connects the in-in formalism with the 
number of vertices in a Feynman diagram. For example, N = always represents propagators at tree level; N = 1 



FIG. 1: The two-point correlation function of primordial tensor fluctuations with one- loop fermion correction involving two 
vertices. 



refers to diagrams with only one vertex; and N = 2 means that there are two vertices involved in each of these 
diagrams; and so on. 

As is well known, in the in-in formalism the Hamiltonian is decomposed into a effectively free part Hq which is of 
quadratic order in fluctuation variables and an interaction part Hj. First, we study the interaction Hamiltonian in 
cubic order. By expanding the action up to next-to-leading order, we get 

S int = I f d A x V=g T» v 5 9llv (25) 



2 . 

Since the scalar and tensor fluctuations are decoupled in tree level, we can write down the perturbed FRW metric 
only involving the tensor modes as follows, 

ds 2 = a 2 (r)(dr 2 - c^dxW) (26) 

where 

9ij = % + % , (27) 

with the tensor perturbation hij satisfying the relations dihij = and ha = 0. To first order in hij, Eq. (|25l) becomes 

S? nt = I /d 4 * a 2 [W(d^) + W(d^) ~ {dii>)lH - (fi^frty] h l3 . (28) 

Similarly, we can expand the action to the forth order and obtain one extra interaction term in order of 0(h 2 ). This 
is a coupling between the mass term of the fermion and two tensor modes, which is given by 

^int = 4 y a3 ' m '4 , ' l Phijhji . (29) 
We will discuss loop corrections arisen from these interaction terms one by one in later context. 



IV. FERMION LOOP CORRECTIONS TO CORRELATION FUNCTION OF TENSOR 

FLUCTUATIONS 

In the previous section, we have derived out two interaction terms (|28|) and (|29p which correspond to the cases of 
two fcrmions coupling to one tensor mode and two tensor modes, respectively. Thus, one can immediately read that 
there exist two different Fcynman diagrams with one-loop correction which are depicted in the following figures. 



A. A Massless Fermion 



As an exercise, we first consider that the fermion is massless. This is because at early times of our universe, the 
energy scale is extremely high and the masses of most particles are much lighter than that of inflaton. In this case, the 
fermion loop corrections can only be contributed by Fig. [T] since the interaction appeared in Eq. (|29[) automatically 
vanishes when m = 0. 



FIG. 2: The two-point correlation function of primordial tensor fluctuations with one-loop fermion correction only involving 
single vertex. 



1. A illustrate of one-loop correction 
We first focus on one interaction term appeared in Eq. ([25]) . which leads to the following Interaction Hamiltonian, 

H i = \f <^ x a 2 {d l ^) 1 ^h l] . (30) 

According to the formula Eq. (|24p . the correlation function with one-loop fermion correction as depicted in Fig. Q] 
can be written as 

(/w(x,t)/i mn (x',t)) 2 

ft rt 2 

dt 2 / dti([fr 7 (ti),[ J Er/(i2),ft m „(t)/i mn (*)]]) 



f 

32 



— oo J — oo 

d 3 Xl d 3 x 2 I dt 2 [ d*ip'V Re 

J — oo J — oo 



x ((h k i(ti)hij(t 2 )h mn (t)h mn (t)) - (hki(ti)h mn (t)h m n(t)hij(t 2 ))) 



(31) 



Note that, in the above formula there is no extra minus sign from the partial operator in Eq. (|30[) . since the contractions 
always take place between the modes e lpx and e _?p x . After contracting the particle modes according to Eq. (|31l) . 
then we get the following three useful expressions: 

(hki (h)hij(t2) h mn (t) h mn (<)) 

= 2 / dWgiE eiqi ' (xl ~ x)+iqHx2 ~ x 'WM 

A 

x h qi (ti)h* qi (t)h q[ (t 2 )h* q ,(t) , (32) 

(hki (*i ) h mn (t) h mn (t)hij(t 2 )) 
= 2 f dVdViE^ 1 '^^' 1 '^ ^^ 

3 A 

x h qi (h)hl(t)h ql (t)h* q>i (t 2 ) , (33) 
= J dVV e ,(p+p ' Hxi - X!) E^4«i)^ s fe)^y(( 2 )^,,(ti). (34) 
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To substitute Eqs. (|32 l 133 1 134 |l into E q- UU, we finally have 

d 3 xe"^ x - x ') (h mn (x, i)/i ro „(x', t)) a 
( 27 3 T T 2 ,_o„ /• J 3« J 3^x3,_ , _/ , „ N p'Vtr( 7 ( M 



f dr 2 / 2 dnRe / d 3 pdV<5 3 (P + p' + q)- 

J — OO J — OO ■/ 



64 J.oo J.oo J pp' 

x ^ e w ($, A)4„(g, A)e mn (g, A'K^g, AOe"^'^^'^ 

A, A' 

x h q {tx)h* q (t) [h q (t 2 )h*(t) - h q (t)h*(t 2 )} , (35) 

where we have made use of conformal time r instead of the cosmic time t via r = f* — . In addition, we have also 
applied a useful relation of the polarization tensor which is given by, 

4'(«> A ) = e y( _ ?> A ) ■ ( 36 ) 
Notice that, one can formulate the Fourier modes of a massless fermion as follows, 

Xp ' s(t) = -^kt) Up ' S ' ^ p ' s( * ) = ^) V ^' ^ 



where 



with u° and v° in the Minkowski space normalized as 



u p ,. = u°e-* T I Wp, s =u°e ipT , (38) 



X! u P:S u p,s ~ X^^s ~~ 2(27r) 3 p 

Since the fluctuations in the interacting picture are viewed as free fields, h q (t) is positive- frequency mode of which 
the form is provided in Eq. (|23p . Integrating over the conformal time in Eq. ((33)) yields the final result: 



d 3 xe l ^ x - x '\h mn ( X ,t)h mn (x',t)) 



2 



■ y d 3 pdV<5 3 (p + p' + q) 



16(27r) 3 M4,/ ^ ^ ^ r ^ y pp' 

tf* -1_ 

<7 6 2< ? (p + p' + < ? ) 



x £ki{q, A) £ ; m (g, A)e mn (g, A'K^g, \')— — — — - . (40) 



A, A' 

With the useful normalization Eq. (|17[) , one can easily do the contraction between momentum and tensor polarization. 

2. The Complete Hamiltonian 

After having studied the one-loop correction from a single interaction term between the massless fermion and 
gravitational waves, then we extend the interaction Hamiltonian into the complete form: 

H I = - l -Jd 3 x aW(c^) + H{dii>) ~ {dii>)lH ~ (d^h^}h tJ , (41) 

which is only contributed by the perturbed action at cubic order. After a process of lengthy computation, we obtain 
the two-point correlation function for tensor modes as follows, 



d 3 xe^~^ (/w(x, i)/w(x', t))a 
1 / dW 3 ( P + p' + q) MVMW - *W - **) 



4(2^) 3 M4 7 " ^ vr ' r ^ W ' 
x J2 e ki(q, A)< Tm (g, A)e Tra »(g, Y)e^-(g, X')— — - - - - , (42) 
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where q is the momentum of the external line in Fig. [TJ and p, p' arc momenta corresponding to the two internal 
lines in the same figure. After contracting all the Lorentz indices, Eq. (|42p can be further simplified as 

J ^xe^--'\h mn {^t)h mn ^,t)) 2 = - 2 Z) 3 ( MV X ^ m ' (43) 



where 



fC(q) = dp dp 

JO J\p-q\ 



2(27r) 3 M4 g 
p+q , (p' 2 - p 2 - q 2 )(p 2 - p 12 - q 2 ) 



\ P -q\ q 4 {p + p' + q) 



A 2 2 ( 12 2 2\2 

4p q — [p — p — q ) 



(44) 



and t q denotes the time scale when the tensor mode with comoving wave number q crosses the Hubble radius. 

By virtue of the dimensional rcgularization analysis, the remaining momentum integral will lead to the following 
correspondence: 

2-7T 

—K{q) -+ q i+s F(S) , (45) 
where 5 denotes the dimensional difference. When 5 — > 0, we have 

JC(q)=^q 5 (F lnq + L) , (46) 
with L being a divergent constant. In order to calculate the coefficient Fq, we differentiate lC(q) and get, 

fl, = ^ . (47, 

In addition, the dimensional difference 8 brings a log correction to the solutions of fcrmion and tensor modes, X PfS 
and h q (see [35| for detailed study). To take into account this effect, and to combine the analysis of dimensional 
rcgularization, we eventually obtain the two-point correlation function of primordial tensor fluctuations under one- 
loop massless fermion corrections, which is given by, 



d 3 X e^~^(h mn ( X ,t)h mn ( X ',t)) 2 H- - ilT H ,}!l\ x ln (— ) . ( 48 ) 



15(2ir) 3 M^q 3 " " V M 

where [i is a physical renormalization scale. Note that, the form of our result is similar to the result of Ref. [24[ in 
which the author studied the loop correction of a fermion field to the primordial curvature perturbation. Our result 
shows that the amplitude of the fermionic loop correction to primordial tensor fluctuation is much larger than that 
to the scalar type curvature perturbation. In addition, we notice that our result is obviously invariant under the 
following rescaling: a — > Xa. x — > x/X, and k — > Xk. This conclusion is consistent with the analysis made in [35| . 

From our result, we notice that the sign of the loop correction arisen from a massless fermion is the same as that 
of a massless scalar field. This conclusion is in agreement with the result of Ref. [H) , in which the loop correction of 
a fermion field to curvature perturbation was analyzed. Here we would like to comment a little bit on this issue. It is 
well known that, there is always a negative sign contributed by a fermionic loop correction in quantum field theory. 
In quantum field theory the loop contribution is calculated based on the in- out formalism and the Wick's theorem 
guarantees that for a closed fermion loop, and thus we have the following relation, 



4>1 ^ltp2 ^2lp3 ■■■ ^n-ltpn 4>n = ~tr 



4>l1p2 i>2i>3 ■■■ Ipn-llpn 4>l4>n 



(49) 



which explains the origin of a negative sign of the fermion loop. However, the perturbation theory in the cosmological 
background is based on the in-in formalism, and correspondingly the contraction relation is replaced by 

^lV'l^2V , 2^3---V'n-l^nV'n = V'l^2^2^3"-V'n-lV'n'0nV'l , (50) 
v v ' 

since there is no time-order in this formalism. 



B. A Massive fermion 



After having analyzed the specific case of massless fermion loop corrections to the correlation function of primordial 
tensor fluctuations, we continue to study the loop effects brought by a massive fermion in this subsection. In the 
case of massive fermion, there exists another interaction due to the expansion of the mass term in the action. This 
interaction could bring a new loop correction with single vertex as shown in Fig. [5] 
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1. One-loop correction with two vertex 



We first consider the same interaction Hamiltonian as what has been studied in the case of a massless fermion. Its 
Feynmann diagram is shown in (|T|). We can simply repeat the process of previous calculation, but need to replace 
the plane wave-function of the massless fermion by the Hankel function which is the mode function for the massive 
fermion field. 

In this case, the correlation function of primordial tensor fluctuations with one-loop fermion correction can be 
expressed as, 



1 J d 3 xe lq ( x " x ') J d 3 x 1 d 3 x 2 J dt 2 a 2 (t 2 ) J * dha 2 l (h)p' k Y Re J d 3 pdVe t(p+p ' ) (xi ~ X2) 
x Y.^ X vAh)X v AWW v ^At 2 )W v ^At 1 ) J d^xd-^^e^'^-^+^-tx.-x') 

s.s' X 

x e fci (gi, A)e^„(qi, A)e mn (-gi,A04(-gi,A')V(*i) /l 9i( < )[V^ 2 ) /l «;^) ~ Vi (k)] > 

(51) 

where the mode function X p s is given by (|22|) . Making use of a little mathematical features of Hankel functions and 
the Dirac algebra, one can get the following useful relation, 



2|r(^)| 4 / ur 27rm 



(27r)8o 3 (ti )a 3 (t 2 ) V 



f C0flh[— ]-lJ 



(P ■ P')(p ™p' + cc.) 



6 jl + 0(e jl ) (52) 



where we have introduced the index [i = | — i@ with m being the mass of the fermion. The second term 0{e > ) in rhs 
of (|52[) is asymmetric and thus will vanish due to symmetry. In addition, we have assumed the perturbation modes 
are super-Hubble scale where the Hankel function is approximately a power law function. Therefore, the following 
calculation is based on the assumption that all the modes of interests are super-Hubble. 

Then we insert the relation (|52")l into (|5"Tj) . After a process of lengthy calculation, we further obtain the simplified 
form of the two-point correlator as follows, 

4ff 4 |r( M )| 4 f_ ul 2nm, \ f q [ p+q A , u 2 2 lj2 2 



ff 4 |r(/i)| 4 / , r 2^™ n \ r q , r p+q , . ri 2 2 , 
(2^^( c ° sh[ TT 1 - 1 )/o d!, L* [4p¥ " C! '' 

2pp' — (p ■ p'Xp^p'" 2 ^ 1 + c.c.) 



r- <? 2 ) 2 1 



(53) 



where we have contracted the indices of tensor modes in detailed calculation. In addition, we can make use of the 
useful relation: |T(/i)| 2 = 7r/cosh[^]. Then we do the integral of ([53)) very carefully, and find the final result is 
composed of some regular functions and many annoying hypergeometric functions. 

In order to grasp the physics of this integral, we would like to express the final result in approximate forms under 
two different limits. First, in the limit of small q value and small fermion mass to, the result of the integral is given 
by 

d 3 xe i(l < X - X '\h m n(x,t)h m n(^,t)) 2 \ q<a H, m< H " ^jl^^j^j^ ' ^ 

up to the leading order. Second, we still consider q to be small but a much larger value of m, and then we get the 
asymptotic form as follow, 

Qui ".inM 27rm l tnnf "" ! 

dW"(»~ \h mn ^t)h mn ^ Mq<<aH , m>>H * ^L}* 1 • (55) 
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From the above result, we can get the following information. For the case of a massive fermion, its one- loop 
correction to the correlation function of primordial tensor fluctuations does not change the momentum dependence 
and thus the primordial power spectrum is always scale-invariant. Moreover, when the fermion mass is very small, 
the contribution of the fermion loop is suppressed by the factor m /H , while the fermion mass is very large, the 
contribution of the fermion loop is suppressed exponentially due to the factor l/cosh 2 [^]. Therefore, one may 
conclude that the logarithmic divergence of the fermion loop appeared in the case of massless fermion might be 
removed by the mass term. However, we should keep in mind that this result is based on the assumption that all 
primordial perturbation modes of interest arc on supcr-Hubble scale and thus we have neglected the contribution of 
the loop integrals inside the Hubble scale during inflation. It is interesting to go to details on this issue and dig out 
the relation between the logarithmic divergence and the loop integrals related to sub-Hubble modes. We would like 
to leave it to future study. 



2. One-loop correction with single vertex 



Then we turn our attention to the calculation of one-loop correction of a massive fermion field with a single vertex. 
The relevant Feynmann diagram is shown in Fig. [2] and its interaction term is given by ((29)) . 

According to the formula developed in Ref. [15|], as there is only one vertex in Fig. @, it corresponds to the N = 1 
term to Eq. (|24p with one communicator in the vacuum expectation. Thus, we have the following expression 

(ft. OT „(x,t)/i ro „(x',i))i 

dii([lfi(ii),/i mn (i)/i mn (t)]) 
t 

dh a 3 (ii) m 

-OO 

(h mn (t)h mn (t)hij(tx)hji(ti)) - (hij(ti)hji(tx)h mn (t)h mn (t)) . (56) 



- / d irid x 2 



Recall that we have (ip(ti)ijj{ti)) = J d 3 pJ2 s W p , s (ti)W PyS (ti) for the fermion field. Again, we make use of the 
Dirac algebra and contract the indices of tensor modes. After a little lengthy computation, we eventually get the 
result of the above integral: 

d 3 xe il(*-~') < h mn {*,t)h mn {x! ,t) > 1 



647rmiJ 4 
(27r) 4 Af4 g f 

64mAfj V 



Im 



/T ra(Ti)A uv 
dT 1 a{T 1 )e~ 2iqTl+2iqr / d 3 p 1 

OO J 



C + ln(— ) 



(57) 



where we have used aH = q at the moment of Hubble-crossing. Ajjy is a physical cutoff at UV regime, and C is some 
integral constant which is of 0(1). 

From the result (|57[) . one can find that the divergence of leading order is the UV divergence, and there is a mixture 
between the UV and IR divergence due to the logarithmic term at the next-to-leading order. Moreover, this loop 
correction is proportional to the factor jj and thus becomes negligible when the fermion field is much lighter than 
the inflaton. Eventually, since the final integral is proportional to q~ 3 , it implies that the scale invariance of the 
primordial tensor fluctuations would not be changed under the fermion loop correction. 



V. CONCLUSION 



To conclude, we in the current paper have studied the loop corrections of a fermion field to the primordial power 
spectrum of tensor fluctuations in the frame of inflationary cosmology. Our calculation is based on the in-in formalism 
and the spinor field is assumed to be a Dirac fermion. In this setup, we have separately considered the loop contribution 
of a massless fermion and that of a massive fermion. For the case of a massless fermion, we obtain a logarithmic 
correction in the correlation function of primordial tensor fluctuations which involves a regularization parameter in 
IR regime. This result is in agreement with previous works fl5l |24||. Notice that, our result shows that this IR issue 
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is irrelevant to the scale factor rescaling. Thus the IR regularization parameter is a physical one, which is consistent 
with the analysis in Ref. (35| . 

Further, we have extended our calculation to the case of a massive fermion field. We find that the inclusion of a 
fermion's mass can be helpful to relax the IR issue since the logarithmic function obtained in the massless fermion 
loop is replaced by a cosh function which is exponentially suppressed at large mass limit. Moreover, at the small mass 
limit, the fermion correction is also suppressed since the loop integral is proportional to m 2 / H 2 . However, we should 
be aware of that this result is based on the assumption that in the loop integral most of the perturbation modes are 
in the supcr-Hubble scale. It is interesting to study in detail the contribution of the sub-Hubblc modes in the loop 
integral and see how it would affect the IR issue in the case of a massive fermion field. We would like to leave it as a 
follow-up project. 

Moreover, due to the fermion's mass term, one can obtain a fermion loop to the correlation function of primordial 
tensor fluctuations with single vertex. In usual quantum field theory, this term can also appear but only contribute to 
the UV renormalization. However, in the context of inflationary cosmology, we find the leading term of the one-loop 
correction is proportional to Afj V while the next-to-leading term is a mixture of the UV and IR divergences. It is 
necessary to continue the study in this issue and see if there is any implication on the relation between inflation 
models and quantum gravity. We notice that the authors of Ref. [3 a ] also acquired a similar result on primordial 
curvature perturbation in the presence of loop corrections of a second entropy field. 

As an ending remark, we would like to point out all the above loop corrections do not change the scale dependence 
of the correlation function of primordial tensor modes. Therefore, we expect a nearly scale-invariant power spectrum 
of primordial tensor fluctuations. 
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